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ABSTRACT 

We report on our second-epoch VLBI and VLA observations of the Type Ib/c supernova 2001em, 
five years after the explosion. It was suggested that SN 2001em might be a jet-driven gamma ray burst 
(GRB), with the jet oriented near the plane of the sky, which would entail relativistic expansion or 
motion. Our VLBI image shows that SN 2001em is still unresolved five years after the explosion. For a 
distance of 83 Mpc (i?o = 70 km s _1 Mpc -1 ), the nominal expansion velocity is 5800±10, 000 km s _1 , 
and the proper motion is 33, 000 ± 34, 000 km s" 1 . Our values are inconsistent with either relativistic 
expansion or motion, but are consistent with the non-relativistic expansion speeds and small proper 
motions seen in other supernovae. In particular these values are consistent with radio emission from 
SN 2001em being due to normal, non-relativistic supernova ejecta interacting with the circumstellar 
medium. Our VLA observations show a power-law decay in flux density since the time of the peak in 
the 8.4 GHz radio lightcurve in -2003. 

Subject headings: supernovae: individual (SN2001em), radio continuum: general, gamma rays: bursts 
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1. INTRODUCTION 

Supern ova 2001em was disco vered on 2001 Septem- 
ber 15 (jPapenkova fc Lil 1200 If ), in the galaxy UGC 
11794 (NGC 7112; z = 0.01935, D = 83 Mpc for 
Hq = 70 km s -1 Mpc" 1 ). It was classified as a 
Type Ib/c on the basis of a n early optical spectrum 
(jFilippenko fc Chornockl[2"001h . Subsequently, however, 
it went on to show strong signs of c ircumstellar in- 
teraction (see IChugai fc Chevalier] 120061 and references 
therein). Among these signs were radio and X-ray lu- 
minosities unprecedented fo r a T ype Ib/c supernova af- 
ter ~ 2 yr (|Stockdale et alj|2004l200l (Poolev fc Lewinl 
12001 . 

It was suggested that the radio emission from SN 
200 lem might be due to the interaction with the cir- 
cumstellar medium (CS M) of a relativistic jet, oriented 
near the plane of the sky (jGranot fc Ramirez-Ruizll2004l . 
see also Paczyhski 2001). Such jets are thought to power 
gamma-ray burst (GRB) events; in particular long dura- 
tion GRB events, which seem to be associated with Type 
Ib/c supernova. 

A prediction of the relativistic jet model for SN 2001em 
is that the source should be resolvable by very long base- 
line interferometry (VLBI) . Along wi th two other groups, 
we u ndertook such observations (jBictcnho lz fc Bartel 
l2005l referred to hereafter as iPaper If) . The two 
other group s reported their res ults in IStockdale et al.l 
(|2005h and iParagi et al.l (|2005h . SN 2001em was at 
best marginally resolved in all three experiments. We 
obtained the strongest limit on the expansion speed, 
namely a 3cr lim it of 70,000 km s _1 on one-sided ex- 
pansion (|Paper J ), which does not support the relativis- 
tic jet model. An altern ate model was proposed by 
IChugai fc Chevalierl (|2006l ). which involves the interac- 
tion of normal (non-relativistic) Type Ib/c ejecta with a 
massive and dense circumstellar shell, produced by mass 
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loss of the progenitor. 

The relativistic jet model, however, is not entirely ex- 
cluded by the small size observed for the radio emis- 
sion region: it is possible that only the "working end" 
of the relativistic jet is radio-bright. In this case, how- 
ever, the radio-bright spot should display a relativistic 
proper motion. At 83 Mpc, a speed of c corresponds to 
0.8 mas yr -1 , which is resolvable by current VLBI arrays 
in less than one year. 

To better measure the expansion speed, and to deter- 
mine the proper motion of SN 2001 we undertook "second 
epoch" VLBI observations in 2006 May, 4.7 yr after the 
explosion, and 1.4 yr after our previous observations. We 
report the results of these observations here. 

2. OBSERVATIONS 



As in lPaper^ . we observed SN 2001em using the "High 
Sensitivity Array," which consists of NRA0 3 's Very Long 
Baseline Array (VLBA; each of 25 m diameter), phased 
Very Large Array, (VLA; 130 m equivalent diameter), 
and Robert C. Byrd telescopes (—105 m diameter), in 
addition to the Arecibo Radio Telescope 4 (305 m diam- 
eter) and the Effelsberg Radio Telescope 5 (100 m diam- 
eter). The VLBI observations were carried out on 2006 
May 27, with a total time of 12 hrs. We observed at a 
frequency of 8.4 GHz with a bandwidth of 32 MHz and 
recorded both senses of circular polarization with a total 
recording bit rate of 256 Mbits s _1 . The data were corre- 
lated with the NRAO VLBA processor at Socorro. The 
analysis was carried out with NRAO's Astronomical Im- 
age Processing System (AIPS). The initial flux density 
calibration was done through measurements of the sys- 
tem temperature at each telescope, and then improved 

3 The National Radio Astronomy Observatory, NRAO, is oper- 
ated under license by Associated Universities, Inc., under cooper- 
ative agreement with National Science Foundation, NSF. 

4 The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which is operated by Cornell University 
under a cooperative agreement with the NSF. 

5 The 100 m telescope at Effelsberg is operated by the Max- 
Planck-Institut fur Radioastronomie in Bonn, Germany. 
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through self-calibration of the reference sources. Finally, 
in addition to using the phased VLA, as an element in 
our VLBI array, we also used the VLA, which was in the 
AB configuration, as a stand-alone interferometer to de- 
termine the integrated flux density of SN 200 lem at 1.5, 
8.4 and 22 GHz. 

We phase-referenced our VL BI obse rvations to the 
same two sources we used in [Paper J . namely JVAS 
J2145+1115 and JVAS J2139+1423. For most of the 
observing session, we phase-referenced to our primary 
phase-reference source, J2145+1115, which is 1.4° away 
on the sky. In addition, we had several ~20 min peri- 
ods where we phase-referenced to our secondary source 
J2139+1423, which is 2.1° away on the sky, and is a 
"defining sourc e" of the Interna tional Celestial Reference 
Frame (ICRF; iMa et al.lfl99l ). The use of the second 
reference source provides both a check on the astrome- 
try and on possible source motions of J2145+1115, and 
allows us to determine coordinates for SN 2001em tied to 
the ICRF. For the observations, we used a cycle-time of 
~240 s, with ~170 s spent on SN 2001em and -70 s spent 
on the calibrator, either J2145+1115 or J2139+1423. 

3. RESULTS 

3.1. VLA Total Flux Densities and the Radio 

Lightcurve 

Using the VLA as a stand-alone interferometer, we 
measured the integrated flux density of SN 200 lem on 
2006 May 27 to be 1.05 mJy at both 1.4 and 8.4 GHz, 
with uncertainties of 0.06 and 0.22 mJy, respectively. At 
22.5 GHz, we did not detect SN 2001em, and instead 
place a 3a upper limit on its flux density of 0.8 mJy. We 
show the radio lightcurve of SN 2001em in Fig. [TJ The 
lightcurve reaches a peak at time after the explosion, t, 
of ~2.8 yr, and shows an exponential decay thereafter, 
with the flux density being proportional to i- 105 ± 02 . 

3.2. Astrometric Results: the Proper Motion 

The astrometric results were obtained by fitting mod- 
els in the u-v plane to the phase-referenced visibility data 
from the 2006 May as well as the 2004 November ob- 
servations. Arecibo was not used for deriving the as- 
trometric results, as systematic phase-shifts between the 
calibrator sources and SN 2001em were observed on base- 
lines involving Arecibo. For the 2004 epoch, we also did 
not use data from the Kitt Peak VLBA station, as they 
were of poor quality. For data from the remainder of the 
array, we treated the sections of data phase-referenced 
to J2145+1115 separately from those phase-referenced 
to J2139+1423. The model-fitting involves the least- 
squares fitting of geometrical models to the VLBI vis- 
ibility data in the Fourier transf orm or u-v plane, and is 
des cribed in mor e detai l in, e.g.. lBietenholz et al.l (|2001f ) 
and lBartel et alj (|2002D . 

We fit a model consisting of a circular Gaussian to 
the fully calibrated visibility data. Since the source is 
only marginally resolved, the choice of model has little 
effect on our astrometric results. For example, fitting ei- 
ther a point-source, an elongated elliptical Gaussian, or a 
spherical shell model produces the same center position 
to well within the uncertainties. The statistical uncer- 
tainties in the fit positions are on the order of 50 //as, 
but the true uncertainties are likely dominated by sys- 
tematic effects, such as those of unmodelled parts of the 



troposphere and ionosphere and evolution of the refer- 
ence sources. We will therefore use an estimated la un- 
certainty of 100 and 280 /xas in each coordinate when 
phase-referencing to J2145+1115 and J2139+1423, re- 
spectively. The uncertainties when phase-referencing to 
J2139+1423 are higher both because it is farther away on 
the sky and because we spent only a fraction of our ob- 
serving time using it as the reference source. For further 
discussion of the limitations on t he astrometric ac c uracy 
obtain a ble with VLBI, see, e.g., iBietenholz et al.l (|2001l . 
I2004blk iPradel et all (l2006h . 

If we use the larger segment of our data in which we 
phase-referenced to J2145+1115 we find a proper motion 
of SN 2001em relative to that source of —66 ± 93 and 
60 ± 93 ^as yr -1 in RA and decl., respectively, or 89 ± 
93 /ias yr~* in total (Fig. [2] shows the locations of the fit 
positions for both epochs on the 2006 image). If, instead, 
we use the smaller portion of our data, in which we phase- 
reference to our second reference source, J2139+1423, 
we obtain a consistent, but less accurate, relative proper 
motion of 192 ± 280 and 20 ± 280 /ias yr" 1 in RA and 
deck, respectively. 

To what extent is our proper motion of SN 2001em 
dependent on the stationarity of the reference sources? 
Since the relative proper motion between SN 2001em and 
J2145+1145 was low, it is likely that the proper mo- 
tion of J2145+1145 is also low, since it would be un- 
likely that both SN 2001em and J2145+1145 have large 
proper motions which fortuitously cancel in our differen- 
tial measurement. It is more likely that the proper mo- 
tion of SN 200 lem is in fact consistent with zero at the 
level of our uncertainties. Our second reference source, 
J2132+1423, is an ICRF source whose proper motion 
has b een measured to be < 40 ^as yr -1 ( Feis sel- Vernierl 
2003). Our measured proper motion of SN 2001em rel- 
ative to it was also consistent with zero to within la. 
We therefore consider it reasonable to assume that our 
primary reference source, J2145+1145, is stationary, and 
that the proper motion of SN 2001em is 

[X = 89 ± 93 ^as yr" 1 , (p.a. = -138°). 

At 83 Mpc, this proper motion corresponds to a non- 
relativistic speed of 

v = 33, 000 ± 34, 000 km s" 1 , (p.a. = -138°). 

We note that even our less accurate value of the proper 
motion of SN 2001em with respect to J2139+1423 is in- 
consistent with a projected speed of c at the 2.2er level. 

3.3. Images and the Size of the Radio- Emitting Region 

To obtain the best imaging results, we phase self- 
calibrated the data for SN 2001em, now including the 
data for Arecibo and combining the sections phase- 
referenced to J2145+1115 and J2139+1423. This self- 
calibration serves to correct the mentioned systematic 
phase-shift seen at Arecibo, and also any phase shifts 
due to possible inconsistencies in the positions of the two 
reference sources. 

In Figure [2] we show the VLBI image of SN 2001em 
on 2006 May 27. No structure is apparent. The source 
is, as it was in our 2004 observations, at best marginally 
resolved at our resolution of 1.7 x 0.8 mas at p.a. — 1° 
(FWHM of a Gaussian; p.a. is east of north). In partic- 
ular, no structure larger than ~ 1.7 mas or ~ 2xl0 18 cm 
is visible. 
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To get a more accurate determination of the size, we 
turned again to model fitting in the u-v plane, which 
has the advantage of giving results independent of the 
convolving beam. Since the source is at best marginally 
resolved, we must choose t he model a priori. We use the 
same models as in lPaperU namely an elliptical Gaussian, 
which we choose to represent a possibly elongated source 
such as a jet, and a spherical shell model appropriate for a 
supernova. The shell model consists of the projection of a 
spherically symmetrical, optically thin shell of emission, 
with the shell width being 20% of the ou ter radius (as was 
found appropriate for SN 1993J; see e.g. lBietenholz et all 
[2003]) . 

We discus sed the u ncertaintie s in the fitting process 
in detail in iPaper 1 (see also iBietenholz et al.l 120011 : 
iBartel et al.ll2002T i7 To summarize, the uncertainties are 
likely dominated by residual inaccuracies in the complex 
telescope gains. We again use the square root of the 
visibility weights in our fits, which results in a more ro- 
bust fit less dominated by a few very sensitive baselines. 
We estimated the uncertainties by allowing the antenna 
gains to vary and by comparing the fitted values with 
and without Arecibo, our most sensitive antenna. 

Using an elliptical Gaussian model appropriate for a 
jet, we obtained a FWHM major axis size of 0.16 ± 
0.18 mas, corresponding to a 3a upper limit of 0.7 mas. 
At 83 Mpc, this limit corresponds to a size of 8x 10 17 cm, 
assuming undecelerated, one-sided, motion since the ex- 
plosion, to a speed of 54,000 km s -1 . For the spherical 
shell model appropriate for a supernova, we find an outer 
angular radius of 0.07 ± 0.12 mas, with a corresponding 
3cr upper limit of 0.42 mas 6 . At 83 Mpc, this corresponds 
to a size of (9 ± 15)xl0 16 cm and an average expansion 
velocity (since the explosion) of 5800 ± 10,000 km s . 

4. DISCUSSION 

We have made a second-epoch VLBI image of 
SN 2001em, one of the first Type Ib/c supernova to be 
observed with VLBI. SN 2001em had an unusually de- 
layed onset of radio emission, and it had been suggested 
that the radio emission might be due to a relativistic 
jet not aligned with th e line of sight, perhaps linkin g 
SN 2001em to a GRB (jGranot fc Ramirez-Ruizl 120041 ) . 
However, an alternate model, where the radio emission 
is due to the normal Type Ib/c ejecta overtaking a dense 
shell left by an episode of strong mass-loss of the pro- 
genitor, can also reproduce the late-onset radio emission 
dChugai fc Chevalierll200l . 

Our present VLBI image (Fig. [2|) shows that 
SN 2001em is still only marginally resolved 4.7 yr af- 
ter the explosion. The radio emission region is not ex- 
panding relativistically. We also measured the proper 
motion of the radio emission region between our 2004 
and 2006 epochs, and found no significant proper mo- 
tion. The 3<7 upper limit on the source's velocity in 
the plane of the sky was 135,000 km s _1 or 0.45 c. The 
radio-bright region, therefore, does not exhibit any rel- 
ativistic projected motion, but rather is consistent with 

6 We note that other models with a geometry different from a 
spherical shell will give different outer radii, but for likely geome- 
tries, the differences to our selected shell model are < 25%. For 
example, a uniform sphere model gives a radius ~ 16% larger, 
whe reas a very thin she ll model gives a radius ~ 10% smaller (see, 
e.g. JBartel et al.|[2002T) . 



being stationary. Different viewing angles can produce 
a wide range of projected speeds even if the true speed 
is near c. We can, however, constrain the range of view- 
ing angles for a standard jet model where the apparent 
jet velocity /3 a pp, the true jet velocity (3 (both in units 
of c), and the angle 9 between the jet and our line of 
sight are related through /3 aDD = j3sa\9 /(I— (3 cos 9) (e.g., 
iPearson fc Zensusl [r987). For example, if (3 — 0.95, then 
only a jet oriented at either 9 > 130° or 9 < 1.4° will 
Papp < 0.45. Thus, if we assume we are viewing the 
Doppler-boosted approaching jet, the only way to rec- 
oncile a true velocity > 0.95 c with our proper motion 
results is to assume that the jet is fortuitously aligned to 
within 1.4° of the line of sight. 

The fact that the radio-emitting region is neither 
expanding nor moving relativistically therefore argues 
against the presence of a relativistic jet in SN 2001em. A 
further argument against such a jet is the detection from 
SN 2001em of strong Ha emissi on with a full w i dth a t 
half maximum of 1800 km s" 1 (jSoderberg et alJl200l . 
Such emission is usually observed in supernovae under- 
going strong circumstellar interaction. 

Both the radio morphology and proper motion 
of SN 2001em, on the other hand, are consis- 
tent with the hypothesis of non-relativistic supernova 
ejecta interacting w i th th e CSM, as was proposed by 
iChugai fc Chevalier! (|2006h . In this model, the late on- 
set radio and X-ray emission is caused by the interaction 
of the supernova ejecta with a dense shell in the CSM 
caused by an episode of rapid mass loss of the progenitor. 
This model can naturally account for the Ha emission, 
which is probably circumstellar gas accelerated by the 
forward shock. 

When we fit a spherical shell model to the radio 
emission, we found a nominal size in 2006.4 of (9 ± 
15)xl0 16 cm and an average expansion speed since the 
explosion of 5800 ± 10,000 km s — 1 . Although our esti- 
mates are rather uncertain, both are compati ble with the 
typic al values usually seen in supernovae (e.g.. lBietenholzl 
l2005h . and in particular with the values used in the 
Chugai & Chevalier model for SN 2001em. We note 
that VLBI observations were also made of Type Ib/c 
SN 2003L, however only a rather large limit of 5.4c 
was obtained for dire ctly measured expansion velocity 
(jSoderberg et al.|[2005f ). In addition, the radio lightcurve 
(Fig. [l} appears to exhibit a power-law decay with flux 
density oc t~ 105 . This is shallower t han that predicted 
for the jet model (jGranot et al.ll2002l ). but such shallow 
decays are seen in some T ype II radio supernova (e.g., 
IBartel fc~l3 ictcnhol3 120031 ) . As an aside, we note that 
our expansion, proper-motion, and flux density measure- 
ments are compatible also wit h a young pulsar nebula, as 
is possibly seen i n SN 1986J (IBietenholz fc Bartell [20061 : 
IBietenholz et alJl200 4a) , although a pulsar nebula would 
not be expected to produce the observed broad Ha emis- 
sion. 

In summary, we obtained second epoch VLBI obser- 
vations of the Type Ib/c supernova 2001em, which had 
unusually late turn-on radio and X-ray emission. We find 
no evidence of relativistic expansion or motion of the ra- 
dio emitting region. Rather, we find that the proper 
motion is consistent with SN 2001em being stationary 
and that the expansion speed is consistent with those of 
normal supernovae. These results favour a model where 
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normal Type Ib/c supernova ejecta are interacting with Research at York University was partly supported by NSERC. 
a dense circumstellar shell. 
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Fig. 1. — The radio lightcurve of SN 2001em at 8.4 GHz, with our additional flux density measurement at 1.4 GHz (shifted slightly for 
clarity) and our 3cr upper limit on the radio flux density at 22 GHz, both for 2006 May. The data points are from VLA measurements 
except that at t = 2.8 yr, which is from the VLBI observations. The error-bars are standard errors consisting of the statistical uncertainty 
and that in the flux density calibration a dded in quadrature. The filled points are our own measurements and the open circles are published 
values from lStockdale et al.l J2005I. [200l) . 
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Fig. 2. — A VLBI image of SN 2001em on 2006 May 27. The peak brightness was 950 ^tjy bm _1 and the rms background noise was 
15 fijy bm — 1 . The contours are drawn at —5, 5, 10, 20, 30, 50, 70 and 90% of the peak brightness, with the lowest contour being at 3c, 
and the 50% contour being emphasized. The FWHM size of the Gaussian restoring beam, which can be compared to the 50% contour of 
the source, was 1.7 X 0.8 mas at p. a. —1°, and is indicated at lower left. The gray scale is labelled in fijy bm -1 . The two crosses show the 
fit position for the two epochs, with the right cross representing the current epoch (2006.4) and the left one the 2004.9 epoch. The origin 
is at R.A. = 21 h 42 m 23?60938, dccl. = 12° 29' 50'.'3001 (J2000). North is up and cast to the left. 



